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Abstract 

Compressed sensing was proposed by E. J. Candes, J. Romberg, 
T. Tao, and D. Donoho for efficient sampling of sparse signals in 2006 
and has vast applications in signal processing. The expicit restricted 
isometry property (RIP) measurement matrices are needed in practice. 
Since 2007 R. DeVore, J. Bourgain et al and R. Calderbank et al have 
given several deterministic cosntrcutions of RIP matrices from various 
mathematical objects. On the other hand the strong coherence prop¬ 
erty of a measurement matrix was introduced by Bajwa and Calder¬ 
bank et al for the recovery of signals under the noisy measuremnt. In 
this paper we propose new explicit construction of real valued RIP 
measurement matrices in compressed sensing from algebraic geometry. 
Our construction indicates that using more general algebraic-geometric 
objects rather than curves (AG codes), RIP measurement matrices in 
compressed sensing can be constructed with much smaller coherence 
and much bigger sparsity orders. The RIP matrices from algebraic 
geometry also have a nice asymptotic bound matching the bound from 
the previous constructions of Bourgain et al and the small-bias sets. 
On the negative side, we prove that the RIP matrices from DeVore’s 
construction, its direct algebraic geometric generalization and one of 
our new construction do not satisfy the strong coherence property. 
However we give a modified version of AG-RIP matrices which satis¬ 
fies the strong coherence property. Therefore the new RIP matrices 
in compressed sensing from our modified algebraic geometric construc¬ 
tion can be used for the recovery of signals from the noisy measurement. 
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University, Hangzhou 310018, Zhejiang Province, China. H.Chen was supported by NSFC 
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Index terms: Compressed sensing, restricted isometry property, 
strong coherence property, algebraic varieties over finite fields, Deligne- 
Lusztig surface, toric surfaces 


1 Introduction 

Compressed sensing is a technique that measures a sparse signal in x G 
which has k non-zero coordinates at sampling rates that are substantially 
much lower than the Nyquits-Shannon rate ( [i EH HU i35i mi Ha IS] ) . For 
a measurement matrix $ which has n budged rows and sampling /c-sparse 
signals in {N columns, n < N), the measurement is y = $ • x, where 
X G R^ is a vector with only at most k non-zero coordinates. We refer 
to |10[ [T^ m] for the effective recovery algorithm of the original signal x 
from y. We say that a matrix $ satishes the restricted isometry property 
(RIP) of order k with the constant 5^ satisfying 0 < < 1, if for ev¬ 

ery /c-sparse vector (that is, only k coordinates of the signal are non-zero) 
X G R^, 1 — ^ ^k- The RIP of the measurement matrix 

$ guarantees the effective recovery of the fe-sparse signal x G R^ from 
y = $ • X (EH EU] via linear programming. Here we should note that the 
RIP is only the sufficient condition of the recovery ( |131ll61ll0] i. It has been 
shown that random matrices satisfy the RIP with high probability. How¬ 
ever in practice the sampling has to be done deterministically. Explicit RIP 
measurement matrices are needed in practical compressed sensing. There 
have been a lot of deterministic constructions of RIP measurement ma¬ 
trices from various mathematical objects such as polynomials over finite 
fields, codes, graphs, and additive combinatorics and number theory( see 
[II El Cl EOl El El El ESI Eg Ea Ell Ea ET]). For further extension of com¬ 
pressed sensing in various practices we refer to la n 0 d ca d- For the 
theoretical analysis of the sparse-order k and the measurement budget n we 
refer to paper m- 

For a, n X N matrix $ with N columns (j)i, , the coherence is 

= maxi^j • We have from the Welch bound. 

It is proved in mm that a matrix $ with the coherence /i# satisfies the 
RIP with the sparsity order k < — 1. Thus it is desirable to give explicit 

construction of matrices with small coherence in compressed sensing. 

The average coherence of a n x matrix $ with N columns </>i,..., </> 7 v £ 
R"’ is defined as Ikj l -HH I analysised in 
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[Si], this parameter is important for the recovery of sparse signals from the 
measurement in the presence of noise. The n x N measurement matrix $ 
satisfies the strong coherence property if and cu# < The 

main result of section 2.2 of [5] gave an algorithm for the recovery of sparse 
signals from noisy measurement provided that the measurement matrix has 
the strong coherence property. 

The first systematic deterministic construction of RIP matrices is due 
to R. DeVore [20]. R. Calderbank and B. Hassibi and their collabora¬ 
tors gave many deterministic constructions of RIP (or statistical isome¬ 
try property) matricse in (S] [29] [28l |9| . In 2011 J. Bourgain, S. Dilworth, 
K.Ford, S. Konyagin and D. Kutzarova [7| gave explicit RIP matrices sat¬ 
isfying n = o(/c^) where n is the row size and k is the sparse order, from 
new estimates about exponential sums and additive combinatorics. Their 
n X N RIP matrices have coherence ^ matching up the asymptotic bound 
= 0^^ niog{n/iogN) )^^^) “ <n< {logN)‘^ (page 149 of |7]). 

From the construction in |6| RIP nx N matrices with size parameters in the 
range {logN)'^-^ < n < {logN)^ and the coherence fi = O (( 
can be constructed. There have been many explicit constructions of RIP ma¬ 
trices from various mathematical objects such as chirp sensing codes m), 
BCH codes mi Reed-Muller codes ( ESI), orthogonal codes (EZj), Reed- 
Solomon codes ( |32j ) and expander graphs ( [29| ). Some early works for other 
motivations [DEI El led to RIP matrices with small coherence. The asymp¬ 
totic behaviour of the coherences of RIP matrices from the constructions in 
[SIE] are nice (see [7], pages 148-149). However the further strong coherence 
property of these RIP matrices has not been considered in these previous 
constructions. 

The basic construction of [20| is as follows. For a polynomial / with 
degree less than or equal to r — 1 in Fp[x] where Fp is a finite field with p 
elements (here p is a prime number, it can also be used for finite field with 
q = elements), the length p^ vector Vf = {f(^a,b)) is determined by its p^ 
coordinates fi^a,b) for (a, 6) G Fp x Fp. Here f(^a,b) = 0 if /(a) / b, f(a,b) = 1 
if /(a) = b. Then the columns of these lenght vectors give a. p^ x p'^ 
matrix. It was proved in [4] that the coherence of this p^ x p^ matrix D 
satisfies that /td < 

This construction has been generalized to an arbitrary projective non¬ 
singular algebraic curve (see ED)- Let X be a projective non-singular al- 
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gebraic curve defined over a finite field Fg of genus 5 , P = {Pi,i^p|} 
be a set of Fg rational points on the curve X and G be a Fg rational di¬ 
visor. The functions are taken from the function space T(G) associated 
with a Fg rational divisor G. For each / E T(G), the length g'|P| vector 
Vf = {f(^a,b)) where (a, 6) E Fg x P, is defined as follows. f(^a,b) is 0 if f{b) / a 
and fafi) is 1 if f{h) = a. Thus the number of the columns is dim{L{G)) 
which can be computed from the Riemann-Roch theorem ([12, 17]). We 
have dim{L{G)) = degG — g + 1 degG > 2g is satisfied. Each function 
/ E T(G) leads to a lenghth g|P| binary vector Vf. This g|P| x dim{L{G)) 
matrix has its coherence fi < (see mm)- 

In this paper we give several constructions of explicit RIP matrices with 
small coherence from general projective algebraic varieties over finite fields. 
The calculation of coherence is based on the counting of rational points of 
these projective algebraic varieties. Thus some results about the estimation 
of the number of rational points of projective algebraic varieties play an 
important role. Then some examples of explicit RIP matrices from Fermat 
surface, projective spaces, ruled surface, Deligne-Lusztig surface and tori 
surfaces are given. The RIP matrices from our construction C can match 
the asymptotic bounds of the RIP matrices constructed in J. Bourgain et 
al m) and | 6 ] in some range. Generally speaking from our new algebraic 
geometric construction RIP measurement matrices with much smaller co¬ 
herence and then much bigger sparsity orders than that of the previously 
constructions can be obtained. On the negative side, we prove that the RIP 
matrices from the DeVore construction and its direct generalization in m 
have no strong coherence property. However we present a modified version of 
our constructsion C (including previous DeVore and AG-code construction 
in m) which leads to RIP matrices satsfying the strong coherence property. 
This shows that explicit RIP matrices from algebraic geometry have small 
coherence and can be used in compressing sensing even in the sampling in 
the presence of noisy. 
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2 RIP matrices with small coherence from general 
projective algebraic varieties 

2.1 Construction A 

Let X be a projective non-singular algebraic curve defined over the finite 
field Fg with q elements, where g is a prime power. Let G = Gi - • -|-Gt be 
a Fg-rational divisor on X and P = Pi • • -|-F]p| be a set of |P| Fg rational 
points on X. For every rational function in the function space L(G), that 
is, every rational function on X with at most order 1 poles at the points 
We associated a length {q + l)(t -|- |P|) vector vj = to /, 
where a G FqU{@},b G {Gi,..., Gt} U {Pi,..., i^P|}. f(^a,b) = 0 if f{b) / a, 
f{a,b) = 1 if f{b) = a and a / @ and /(a,b) = —1 if f{b) = a = @ and b is 
a pole point of the function /. Then we have a (g -|- l)(t -|- |P|) x q^ matrix 
^G.P where I = dim{L{G)) is the dimension of the function space which 
can be computed from the Riemann-Roch theorem. 

Theorem 2.1. The coherence of this matrix p < (^eg(G)+ | P |' 

Proof. For any function / G L(G), the length (g-|-l)(t-|-|P|) vector Vf is 
non-zero only at the points (/(Pi), Pi),..., (/(F)p|),F)p|) and (@,Gii,..., (@,GiiJ 
where Gi^, ...,Gi^^ are the poles of this function /. Thus | < Vf,Vg > \ < 
degG + min{tf,tg}. Here tj and tg are the numbers of pole points of func¬ 
tions / and g. On the other hand ||uj|| = y^|P| +tf and ||vg|| = ^jPl -|- tg. 

From an easy computation in calculus the coherence attains its up¬ 

per bound only when / and g have 1st pole at all t points Gi,..., Gt. Then 
tf = tg = t = degG and \ < Vf,Vg > \ = degG + t = 2degG. In this case 
llwll = ll'^sll = \/|P| + degG. 

We can also get another real valued RIP measurement matrices as fol¬ 
lows. Let X be a projective non-singular algebraic curve defined over the 
finite field Fg. Let G = tG be a Fg-rational divisor on X where G is a Fg 
rational point of X and t is positive integer. Let P = Pi i-]p| be a 

set of |P| rational points on X. For every rational function in the function 
space T(G), that is, every rational function on X with at most order t pole 
at the point G. We associated a length {q + 1)(1 -|- |P|) vector vj = f(^a,b) 
to /, where a G F^ U {@},b G {G} U {Pi,..., P|p|}. f(^a,b) = 0 if f{b) / a, 
f{a,b) = 1 if /(&) = ® b ^ G, f(^a,b) = if 6 = G and / has a order 
m pole at the point G. Then we have a {q + 1)(1 -|- |P|) x g* matrix $g,p- 
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Here I = dim{L{G)) is the dimension of the function space which can be 
computed from the Riemann-Roch theorem. 

Theorem 2.2. The coherence of this matrix p < ' 

Proof. For any function / G R(G), the length (g + 1) x (1 + |P|) vector 
Vf is non-zero only at the points (/(Pi), Pi),(/(P|p|), P|p|) and (@,G). 
Thus \ < Vf,Vg > \ < degG + tftg where tf and tg are the pole orders of the 
functions / and g at the point G . On the other hand ||uj|| = y|P| + 

and ||vg|| = y^|P| -|- tg. From an easy computation in calculus the coherence 

Ih/HHl upper bound only when / and g have t order pole at 

the point G. Then tf = tg = t = degG and | < Vf,Vg > \ = degG + In 
this case ||u/|| = ||ug|| = \/|P| -|- degG'^. 

Remark 2.1 In Construction A, though the coherence is slightly worse 
than that in [201 [31] , it gives us real valued measurement matrix, while only 
binary measurement matrices were given in |2Ul IHT] . 


2.2 Construction B 

Let Pp^ be the projective plane over the finite field Fg, be the set 
of all non-singular plane curves in Pp^ defined over Fg with their degrees 
equal to r. Then is an open dense algebraic subset in the dimension 
^(r -|- l)(r -|- 2) — 1 = ^r(r -|- 3) projective space (Bertini’s theorem, page 
179, [12]). We denote the number of Fg rational points of this set as |Br|. 
For any non-singular plane curve f{x,y,z) = 0 defined over Fg, the length 
q‘^ + q + I vector v(f)^, where h takes over all Fg rational points of the 
projective plane, is defined as follows. v(f)^ is zero if h is not on the curve 
f{x,y,z) = 0 and v(f)^ is 1 if /i is on the curve f{x,y,z) = 0. Thus 
the Euclid norm of this vector is exactly the number of the Fg points of 
this curve f{x,y,z) = 0. Then ||v(f)^|| > q + I — 2gy/q ([7, 17]) , where 
g = ^(r — l)(r — 2) is the genus of this curve. On the other hand from 
Bezout’s theorem, the inner product | < v(fi)^, v(f 2 );j > | < We have a 
(g^ + O' + 1) X |Br| binary matrix The coherence of this binary matrix 

- g+i-(r-l)(r-2)^- 

If we take r = 2 we have the following result. 
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Proposition 2.1. For any prime power q, we have a (( 7 ^ + g + l) x — 
2q^) binary matrix $2 with the coherence 

Proof. We need to prove IB 2 I > q^ — q'^ — 2q^. It is well-known that a 
plane quadrics is singular if and only if the corresponding matrix is not full 
rank. In the case the matrix of a plane quadrics is a symmetric 3x3 matrix. 
These singular plane quadrics are on a degree 3 hypersurface in Pp^. From 
the Segre-Serre-Sorensen bound (Theorem 4.4 [30] and [2l|) , B 2 is the set 
of Pp^ minus a degree 3 hypersurface, this degree 3 hypersurface have at 
most 3g^ + q^ + q'^ + q + 1 Fg rational points. 

If we take r = 3, we need a lower bound for IB 3 I. 

Proposition 2.2. For any prime power q, we have a + q + l) x IB 3 I 
binary matrix $3 with the coherence where IB 3 I > q^ — 6q^. 

Proof. It is well-known that a cubic plane curve is singular if and only 
if its discriminant, which is a degree at most 7 polynomial of coefficients, is 
zero (see page 48-51 of [3l| and appendiiex A). From Segre-Serre-Sorensen 
bound (Theorem 4.4 of [30|), IB 3 I > — 6 g®. 


More generally from the classical resultant theory ( [24j ) it is known that 

the reducible or singular plane curves of degree at most r in Pp ^ are 

ir(r-|-3) _ ^^4 _ 


in a degree at most hypersurface. Thus we have |Br.| > ga 


l)g 


ir(r-|-3) —1 


Proposition 2.3. IFe/ia?;e a (g^-|-g-|-l) x (g 2 ''('’+^) — (r^ —l)g2^F+3) 

2 

matrix with its coherence at most ■ 


This gives us a lot of matrices with small coherences when r is relatively 
small compared to g. However this construction is worse than DeVore con¬ 
struction [ 20 ] . 


We generalize the above construction to the case that Pp is replaced 
by a non-singular algebraic projective surface X defined over Fg. The set 
of all Fq rational points of this surface is denoted by X(Fq). For a very 
ample divisor D on X, we need to use the linear system Lmear(D) which 
consists of all divisors on X linearly equivalent to D ([27], Chapter V). We 
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denote B(D) the algebraic set of all divisors in the linear system Linear{D) 
which are non-singular curves. This is an open dense algebraic subset in 
the projective space of the dimension dim(Linear(D)) (Bertini’s theorem, 
page 179, EZj). In many cases this dimension can be computed from the 
Riemann-Roch theorem of algebraic surfaces m, page 362). From the 
adjunction formula the genus g of these non-singular curves in this linear 
system Linear(D) satishes 2g — 2 = D.(D -|- Kx). The number of intersec¬ 
tion of two such non-singular curves is at most D.D. 

For any non-singular curve / in the linear system Linear(D), the length 
|X(Fq)| vector v(f)^, where h takes over all Fg rational points on the surface 
X, is defined as follows. v(f)^ is zero if h is not on the curve / and v(f)^ is 1 
if h is on the curve /. Thus the Euclid norm of this vector is exactly the num¬ 
ber of the Fg points of this curve /. Thus we have ||v(f)^ 11 > g +1 - 
where 2g — 2 = D.(D -|- Kx) is the genus of these curves ( |27] page 361). 

Theorem 2.3. We have a |X(Fq)| x |B(D)| matrix with its coherence 
at most where g = ^(D.(D + Kx)) + 1. 

In Theorem 2.3 X is taken as the Fermat surface in Pp ^ defined by 

= 0 ([Ml). Then |X(Fq 2 )| = + l){q^ + 1 ) 

m)- Simply we take D as the hyperplane divisor, D.D = {q + 1). How¬ 
ever we know that all hyperplane sections of the Fermat surface have at 
least {q — l)‘^{q -|- 1) rational points (see Lemma 3.1 in section 3), and B(D) 
can be all these hyperplanes. Thus |B(D)| is q^ + q'^ + q'^ + 1. We have 

a {q^ -|- l){q^ -|- 1) x {q^ + q^ + q^ + 1) matrix with the coherence at most 

_ Q+l _ _ 1 

(g-i)^(9+i) ~ {q-W ■ 


2.3 Construction C 

Let Y be a non-singular algebraic projective manifold defined over F^. The 
set of all Fg rational points of this manifold is denoted by Y(Fq). For 
an effective divisor D on Y, we will use the function space L(D) which 
consists of all rational functions on Y with poles at most — D f|27]L In 
many cases the dimension of this function space can be computed from the 
Riemann-Roch theorem ([27]). For any rational function / E L(D), the 
length q ■ |Y(Fq) — D| vector v(f)^, where h = (a, 6 ), b E Y(Fq) — D and 
a E Fg, is defined as follows. v(f)^ is zero if h = (a, h) satisfy f{h) ^ a, and 
v(f)^ is 1 if /i = (a, b) satisfy f{b) = a. Thus the Euclid norm of this vector 
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is exactly |Y(Fq) — D|. The inner product satisfies | < v(fi)^, v(f 2 );j > | is 
at most the number of zero points in Y(Fq) of the function /i — / 2 . That 
is, the absolute value of this inner product is smaller or equal to the maxi¬ 
mal possible number of Fg rational points of members of the linear system 
Linear(D). We denote this number by N(D). 

Theorem 2.4. We have a q ■ |Y(Fq) — D| x qd-im{L{T>)) with the 

coherence at most 

For example Y = Pp^ and D = rH, where H is the hyperplane di¬ 
visor of the projective space. Then |Y(Fq) — D| = q^, dimL(rH) = 

, Y(D) < rg” ^ + q"^ (Segre-Serre-Sorensen 

. ("r] 

bound, ISO]). We have x q\ ' matrix with the coherence at most 
^^ "When n and r are quite close to q this example is 
better than the outputs of the DeVore construction. 



3 Examples 

3.1. Fermat surface. In this example we need the following Lemma. 

Lemma 3.1. Let t < q + 1 be a positive integer satisfying gcd{q^ — 
l,t) = 1. For any degree t hypersurface Yt in Pp., there are at least 

(g — l)^(g -|- 1) Fg 2 rational points in the intersection ofYt and the Fermat 
surface = 0. 

Proof. For any degree t hypersurface defined by a homogeneous poly¬ 
nomial /(Yq, Yi, Y 2 ,X 3 ) = 0 , when we fix (g — l)(g -|- 1 ) possibilities of 
Yq and Yi as and i = 0 , ...,g (here 6 * is a primitive 

element of the multiplicative group F* 2 .), we note that we get two equations 

Y|+^ + Y|+^ = a and atY| + Y3 + • • • + aiY2Y‘-^ + aoX^ = d, 

where a E Fg and d E Fg2. Then at{^Y H-+ ^0 = When Y3 

is changed to Y 30 *('?“^\ i = 0, ...,g, the first equation from Fermat surface 
is satisfied. There are g-|- 1 possibilities of since gcdiq^ — l,t) = 1. Thus 

2 

if f < g -|-1 there are at least ^ • (g -|- 1 ) images and there are at least one so- 


9 




lution for the second equation. On the other hand from the homogeneousity 
we have at least = {q — 1)^(9 + 1) such solutions. 


From Construction B, if the conditions t < q and gcd{t,q^ — 1) = 1 are 
satished we have + l)(g^ + 1) x with the 

* when D = tH is taken, here H is the hy- 


coherence = 




perplane section of the Fermat surface ([M])- Here we should note that the 
linear system Linear{til) on X is the same as the linear system of degree 
t homogeneous polynomials on Pp (ESI. When c^/q <t<q +1, c is a 

H 

suitable positive constant, and gcd{t,q‘^ — 1) = 1, the matrices from this 
construction have smaller coherences than that of matrices from the DeVore 
construction 


3.2. Projective spaces. From construction C, we take Pp then we 

r 


matrix with the coherence at most ^ + 




< 


have p^'’’ X p 

^ ^ when p tends to the infinity. On the other hand if we want to get this 

matrix from Hermitian curve using the construction in [31] (page 5040 of 
/O + r\ 


I 1 p'^—p^ 

■+5-1, 9 = 


The coherence of this matrix from 


[31]) , set s = 

the construction is When c is a small positive real constant satisfying 


s = 



jT ■ 


+ g — 1 < p^. We have 


9 + r 


r®. Set r = cp where c is a 


positive constant satisfying < c < p~^/^ (then p^ < s < p®, Hermitian 

curve RIP matrices have smaller coherences than DeVore construction RIP 
matrices). The ratio f^He.rmiUan _ ^ • When p is a prime number 

and tends to the infinity, this ratio tends to the infinity. This illustrate that 
our construction C is much better than the direct curve-based generalization 
m of the DeVore construction in |20j . 


3.3. Ruled Surface. We take X = Pk x Ph in construction C. 

t q t q 

The counting of rational points of members in a linear system LineariY)) 
on a surface is treated as follows in [26]. If the Fg rational points of X are 
distributed on several Fg rational curves Ci, ...., in X. Then we count 
how many of curves can be in a member of this linear system and count the 
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intersection numbers of the divisor D with these curves Ci, i = 1, In 
this case the set of rational points on Pp^ x Pp_^ is naturally the disjoint 
union of {q + 1) sets of Fg rational points on curves pi x Pp^, where pi, 
i = 1,g + 1 are (g + 1) rational points of Pp^. Thus if we take the divisor 
D of type {di,d 2 ), that is, polynomials f{x,y,z,'w) which are homogeneous 
in x,y with degree di and is homogeneous in z,w with degree d 2 , we get 
a linear system with dimension {di + l)(d 2 + !)• If di + ^2 < g + I, there 
are at most —did 2 + (ii(g + 1) + ^ 2(9 + 1) rational points on any member of 
this linear system. 1[26]1. We have a g^ x g(‘^i+i)('^ 2 +i) j^jp matrix with the 
coherence at most _ 


Set q = p ^ , we have a p^^ xp 3 ('^i+i)('^ 2 +i) jpiatrix with the coherence at 
most ! + 1 ) _ gg|. range < di,d 2 < p^ 

(and then did 2 ~ s < p® as required in [3T]. If di = d 2 = cpa for suitable 

prime number p and suitable positive number c satisfying 0 < c < ps. For 

the RIP p^^ X p('^i+i)('^ 2 +i) rnatrix from Construction C. the coherence is 
10 

{di+d2){p'3-+l)-did2 ^ _2c_1 


If we want to get this matrix from the construction [31] (page 5040) from 

Hermitian curve, the degree of the divisor s = 5 ('^i+b('^ 2 +i) _|_ -^ 

g ^ —I" ^ 2 ^ ■ Then the upper bound of the coherence from Hermitian 
curve construction in m is ^ ~ gp /3 + ^ Then the ratio of two 

coherence is « 5/12c. If < c < Then the coherence of 

/-^Hermitian 0 

our construction C is much better. When p tends to the infinity, the ratio 
tends to 0 . 


3.4. Deligne-Lusztig surface. Deligen-Lusztig varieties were in¬ 
troduced in [18] and had been found very useful in coding theory m) 
since there are many rational points on them. In this section we use the 
Deligne-Lusztig surface X of the type defined over the finite field Fg 2 
as in |55| section 4 pages 542-545 to construct RIP matrices. There are 
(g^ -|- l)(g^ -|- l)(g^ + 1) rational points on X. The divisor L as in [26] section 
4 is used. It is the pull-back of the tH where H is the hyperplane section 
on Pp and t satisfying t < g^ -|- 1 as required in section 4 of [26|. Then 

in the linear system Linear{lAj there are at most t((g^ -|- l)(g^ -|- 1) -|- (g + 
l)(g^ + 1)(9^ + 1 ~ ^)) rational points. The dimension of this linear system 
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cqt^ where c is a positive constant. From 


construction C we have q^{q^ + 1)(q''^ + + 1) x q'^*' RIP matrix with 

the coherence +^) ~ Set g when t 

is in the range < t < we have ^^construction.c tends to 0 when p 

l-^H ermit ian 

tends to the infinity {pHermiUan as in page 5040 of [SI])- 


3.5. Toric surfaces. Toric varieties are typtical geometric objects in 
algebraic geometry and complex differential geometry (1211 and played an 
interesting role in algebraic-geometric coding theory ([IS]). In this section 
we give some RIP matrices from toric surfacses and show that their coher¬ 
ence are quite small. 

Let C be the set of all integral points . We denote 6 a primitive 
element of the finite field Fg. For any integral point m = {mi, m 2 ) € 7? we 
have a fnnction e(m) : F* x F* —>• Fg defined as e(m)(0*, 0^) = fgj. 

z = 0,1,..., q — l and j = 0,1,..., q — 1. Let A C be a convex polyhedron 
with vertices in and L(A) be the function space over Fg spanned by these 
fnnctions e(m) where m takes over all integral points in A. In the following 
cases of convex polyhedrons these functions are linearly independent from 
the result in [25] . 

For each function / G L(A) we have a length q x {q — 1 )^ vector 
v(/) = (/(a,fe)) where (a,6) G F* x F* x Fg defined as follows. f(^a,b) = 0 if 
/(a) ^ b and f(^a,b) = 1 if /(a) = b. Then we have such vectors 

and a q{q — 1 )^ x qd-^m{L(A)) uiatrix The following cases as in the main 
results Theorem 1, 2, 3 of [25] are considered. 

1) A is the convex polytope with the vertices (0,0), (d, 0), (0,d) where d is 
a positive integer satisfying d < q — 1] 

2) A is the convex polytope with the vertices (0,0), (d, 0), (d, e -|- rd), (0, e) 
where d,r,e are positive integers satisfying d < q — 1, e < q — I and 
e + rd < q — 1] 

3) A is the convex polytope with the vertices (0,0), (d, 0), (0, 2d) where d is 
a positive integer satisfying 2d < q — I] 


We have the following result from the main results Theorem 1, 2, 3 of [25] . 


Proposition 3.1. In the above cases the matrix is a RIP matrix 
whose coherence satisfying the following 
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1) ^($a) < in the case 1); 

2) /i($A) < in the case 2); 

3) ^($a) < in the case 4); 

The size parameter dim{L(A)) is the number of the integral points in A. 
That is, 

1) dim{L{A)) = in case 1); 

2) dim{L{A)) = (cZ + l)(e + 1) + in case 2); 

3) dim{L{A)) = + 2(i + 1 in case 3). 

Comparing with Exampel 3.3 for these RIP matrices from toric surfacses 
the coherence is smaller in some range of paramters. Moreover the coherence 
is better than that of RIP matrices from DeVore construction [20] and its 
direct generalization m in some range of parameters. 


4 zhl-Randomized RIP matrices from algebraic ge¬ 
ometry 

We observe that if the I’s coordinates in DeVore and our construction are 
changed to —I randomly, the conclusion about the coherence is still true. 

Proposition 4.1. We get RIP matrices with ^randomized ^ ^ in DeVore 

construction and fJ,randomized < | y(^^^-d | Construction C if the coordi¬ 
nates 1 ’s are changed to ±1 ’s randomly. 

Proof. For any two columns v(f) and v(g) in the constructed matrix, 
the number of the common non-zero coordinates is at most r in DeVore 
construction and V(D) in our construction C. The Euclid norms of these 
columns are not changed. The conclusion follows directly. 

The following result is about the mathematical expectation of coherence 
after this ±I-randomization of DeVore and our C constructions. 

Theorem 4.1. The probability that the coherence of the Al-randomized 
DeVore construction satisfying ^randomized < ^ is at least and the math- 
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ematical expectation of the coherence is 




E(^: 


randomized 


)< 



. Similarly the probability that the coherence of the zizl-randomized Con¬ 
struction C satisfying tTrandomized < | Y(F^)-D | least and the 

mathematical expectation of the coherence is 


7V(D)-4E3r «'• 


E(/irandomized) ^ 



|Y(Fq)-D| 


Proof. Suppose there are L common non-zero coordinates for two 
columns of the constructed matrix v(f) and v(g). Then the possibility 
that the two columns have the same 1 or —1 at one common position is 
and the possibility that the two columns have the different 1 and —1 at this 
position is Thus the possibility that the inner product < v(f),v(g) >= 

±(L —2rc) is f ^ (“)^- The mathematical expectation of the inner product 


E(< v(f),v(g) >) = = L - when 




L is odd and L — ^ 

elusion follows directly. 




I (—)^ when L is even. The con- 

wr2^ 


In random construction of RIP matrices, we only can get RIP matrices 
with high probabilty. In our ±l-randomized RIP matrices from algebraic 
geometry, the constructed matrices satisfy the RIP property deterministi¬ 
cally and the upper bound of coherence can be derived from counting of 
rational points. Moreover if ±1 coordinates are randomized and the sizes 
and sparsity orders are large, the coherences are much smaller in the sense 
of average. Thus in practical signal processing it is more suitable to use 
±l-randomized RIP matrices from algebraic geometry. 


5 Strong coherence property 

In this section we first prove the following two results. 
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Theorem 5.1. The RIP matrices from DeVore construction and its di¬ 
rect AG generalization do not satisfy the strong coherence property. 

Proof. We need to calculate the sum '^deg{f)<r-n^;^- Since ||uj|p = p 
for each polynomial / G Fp[x], '^degU)<r-it^\ = ^'^deg{f)<r-iVf- On the 
other hand for each point (a, 6) G Fp x Fp, there are exactly p^~^ polyno¬ 
mials /’s in Fp[a;] with degrees deg{f) < r — 1 satisfying /(a) = b. We have 
Tideg(f)<r-iVf = where 1 is a length p^ vector with all coordinates 

1. Thus ■ In the DeVore construction The conclusion 

follows from the condition r < p. 

For the direct algebraic-geometric generalization in m, similarly we 
have ||f/|P = |P|. For each {a, Pi), the coordinate of S/eL(G)j]i^ 

{a, Pi) is jf functions in L{G) are zero at Pi. Thus 

^/6i{G)i]^ = Then w# = conclusion 

follows directly from the fact dimL{G) = degG — g 1. 

Theorem 5.2. The constructed measurement matrices in Construction 
C does not satisfy the strong coherence property. 

Proof. Since ||v(f)|p = |Y(Fq)| — D|. On the other hand for each point 
(a, b) in Fg x (Y(Fq) — D), there are exactly qdim(L{T>)-i) rational functions 

/ satisfying j{b) = a. We have ^ ^dimL{u)_i • Thus it has no strong 

coherence property. 

In the construction of DeVore [20] and our construction C, the coordi¬ 
nates of 1 can be changed to ±1 randomly and the conclusion about the 
coherence is valid. We give some results of this ±l-randomized AG con¬ 
struction in the previous section. If the ±1 could be arranged properly, 
this would give us the RIP matrices with small coherence and satisfying the 
strong coherence property. 

Modified AG-RIP matrices in compressed sensing. 

An even number of ±l’s are said balanced if the numbers of the 1 and 
— 1 are equal and an odd number of ±l’s are said balanced if the difference 
of the numbers of 1 and —1 is ±1. 


15 









A. Balanced ztl in each v(f). 

We divide the set Y(Fq) — D to two balanced parts, one part contains 

IY('p' _j)| 

—L_ai-1] points, we color these points as red. The another part contains 

lYfF ) Dl 

I Y(Fq) - D| - points, we color these points as blue. 

B. Balanced ±1 at each point (a, 6) G Fg x (Y(Fq) — D). 

We fix a base of the function space T(D) as /i,..., /t where T = dim{L(D)). 
For each function / G T(D), set / = oi/i + • • • + axfT- Suppose q = p^. 

If p is an odd prime, let Tr : F^ —)■ Fp = {0, ...,p — 1} be the trace func¬ 
tion. The parity is determined depending on the parity of the element in 
{0,1, ...,p — 1}. We will use the parity of the element Tr{ai or) to 

balance the ±l’s at each point (a, 6). 

In the case p = 2 we use the parity of the Hamming weight ret ((Tr (02 ),..., Tr (ot )), 
that is, the number of I’s in (Tr(a 2 ),..., Tr(a'r)). Here without loss of gen¬ 
erality we assume fi{b) / 0. 

Assignment of ±1 for each v(f). 

Suppose p is an odd prime. For each column v(f) and each (a, b) G 
Fq X (Y(Fq) - D) satisfying f{b) = a, we set (-l)^+P<^rity{Tr{ai+-+aT) if 
p is odd and (^—i'j^+po-rity{wt(Tr(a 2 ),...,Tr{aT))) ass um e that /i(6) / 0) if 
p = 2. Here A = 1 if 6 is red and A = —1 if 6 is blue. v(f) = 0 at the point 
(a, 6) G Fq X (Y(Fq) — D) if f{b) / a. Then in each column v(f), it is 
obvious ±l’s are balanced from A. On the other hand we have the following 
Lemma. 

Lemma 5.1. For each point {a,b) G Fq x (Y(Fq) — D) and ra¬ 
tional functions f G T(D) satisfying f{b) = a, the coordinates of 
columns at {a,b) is 1, the ±l’s of other columns are 

balanced. 

Proof. Suppose a5/i(6) a^fT{b) = a, then the set of all vectors 

(oi — a^,..., aT—a^) satisfying ai/i(6)-|-- • ■ Farfrib) = a is a dimension T—1 
subspace F^“^ of T(D) = F^. Therefore the set of all (oi, ■.■,aT) G F^“^ 
satisfying Tr(ai —a^-l-- • ■+aT — a^) = 0 is a dimension s(T—1) —1 space over 
Fp. There are columns such that the coordinates of these columns 

at (a, b) is 1. In the case p is an odd prime, set : Fp = {0, ...,p — !}—)• Fp, 


16 



4>{x) = —X. It is clear that (p changes the parity. In this case the 2nd conclu¬ 
sion follows from the linearity of Tr over Fp. When p = 2, since fi{b) 7 ^ 0, 
the variables 02 , nr are free and the conclusion follows directly. 


Theorem 5.3. The RIP matrices from modified DeVore construction 
satisfy the strong coherence property if r{r — 1) < If cl) dI(D) > 

— and b) dim{L{'D)) < are satisfied, the RIP matrices 

from the modified construction C has the strong coherence property. When 
p is fixed and q = p^ is sufficiently large, this condition a) is valid. 

Proof. From Lemma 5.1 we have — l. 

Thus the average coherence is at most ^ ps 7 ' _i~^ < ^ ■ In the 

modified DeVore construction it is clear that the coherence is exactly 
The 1st conclusion follows directly. In the case of construction C, it is direct 
from our assignment of ±l’s that the coherence is exactly | y^^^d | • 

2nd conclusion follows from Weil-Deligne theorem |19] or Lachaud bound 
[23] directly. 

In the case Y is projective space and D = rH where H is the hyper¬ 
plane section. Then the condition a) in Theorem 5.3 is iV(rH) > ~ . 

This condition is satisfied automatically. The condition b) in Theorem 5.3 is 
g" 

< - - -. This condition can be satisfied when r/q is relatively 

IdOlogq 

small. In Example 3.3 the condition b) in Theorem 5.3 is {di + l)(d 2 + 1) < 
iQQiogq ■ The condition a) is {di -|- d 2 ){q -|- 1) — did 2 > The condition 
a) and b) are satisfied if di and ^2 are 0(g*) where f < 1. In Example 3.4 
Deligne-Lusztig surface case if f < the condition b) is satisfied and the 
condition a) is satisfied automatically. For RIP matrices from toric surfaces 
in 3.5 the condition a) is satisfied automatically and the condition b) is not 
difficult to be satisfied. 

From above examples we can see that the conditions a) and b) in Theo¬ 
rem 5.3 are not difficult to be satished for most algebraic varieties and very 
ample divisors. Thus from our modified algebraic geometric construction 
of RIP matrices, there are many candidates satisfying the small coherence 
property and strong coherence property for the practical use. 


C n + r\ 
r ) 
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6 Asymptotic bound 


Theorem 6.1. From our construction C, we can construct n x N RIP 
matrices from algebraic curves over ¥^2 satisfying the aysmptotic bound 

parameter range logN <n< {logNy. 


Proof. From the Drinfeld-Vladut bound f|36jb there is a family of pro¬ 
jective curves X/i with genus gh and Nn rational points over the fixed finite 
field Fg 2 satisfying = Q- On each such curve Xh we take a rational 

divisor D/j of degree t^gh satisfying 2 < th < q- We consider the w dimen¬ 
sional algebraic variety X/i x • • • x X/^ (w copy product) and the divisor 
D = SX/i X • • • X Dfi X • • • X X/j. Then dzm(L(D)) > {{th — l)gh)^- On the 
other hand there are at most t(;t/i|Xh(Fq 2 rational points on each mem¬ 
ber of the linear system LinearfD). Thus we have a x 

RIP matrices with the coherence at most — '‘iv ^ i • It is 

|-X.hRq2)r |-X-hRq2| 

clear when 2{th — the RIP matrix size is in the 

range logN < n < {logN^. 


On the other hand ~ (gi+fg2"i/3 ~ when gt 

tends to the infinity. | Xh(F^ 2 ) | ~ ^ when gh tends to the infinity. Thus 
when w,th,q are constants the conclusion follows directly. 


7 Summary 

Explicit RIP measurement matrices are needed in practical application of 
compressed sensing to signal processing. In this paper general method of 
contracting explicit RIP matrices from general algebraic varieties over finite 
fields are presented. Many examples of RIP matrices with better perfor- 
mace than previous works are given. We also indicate ±l-randamization of 
these RIP matrices from algebraic geometry are more suitable for practical 
compressed senseing in the presence of noisy. 
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